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The non-equilibrium dynamics of a rapidly quenched spin-1 Bose gas with spin-orbit coupling is 
studied. By solving the stochastic projected Gross-Pitaevskii equation, we show that crystallization 
of merons can occur in a spinor condensate of 87 Rb. Analytic form and stability of the crystal 
structure are given. Likewise, inverted merons can be created in a spin-polarized spinor condensate 
of 23 Na. Our studies provide a chance to explore the fundamental properties of meron-like matter. 
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Spin-orbit coupling (SOC) [1 is an ubiquitous quan- 
tum phenomenon, which links the internal (spin) and or- 
bital (linear or angular momentum) degrees of freedom 
of a particle. The best known example of SOC arises in 
the motion of electrons in the atom, where the electron's 
orbiting around the nucleus can affect the orientation of 
electron's spin. Recently, SOC with tunable strength has 
been realized by Lin et al. [2 by shinning two orthogo- 
nal laser beams intersecting in a pseudo spin- 1/2 Bose- 
Einstein condensate (BEC) of 87 Rb [3]. The lasers are de- 
tuned from Raman resonance so that the momentum and 
spin can couple via exchanging photons in the two orthog- 
onal beams. This opens new possibilities to simulate the 
role of SOC for a wide range of phenomena in condensed 
matter physics by using ultracold atoms, promising ap- 
plications to quantum computing [4], spintronics devices 
[5j [6] and topological insulators [7H2] • Inspired by the ex- 
perimental realization in Ref. [2], theoretical extensions 
using atomic BECs have been studied by a number of au- 
thors, including the direct incorporation of SOC into the 
spin- 1/2, 1 and 2 BECs [T0HT3] , and non-trivial ground 
state structures have been predicted. Novel excitations 
have also been predicted in fermionic gases [H] , such as 
the Rashba pairing bound state (Rashbon) [T5j [16] • 

So far, all studies on the ultracold atoms with SOC 
were focused in the cases of zero temperature. It is fun- 
damentally important to see how the nonlocal nature of 
SOC affects the pattern of spontaneous symmetry break- 
ing. In this work, we study the non-equilibrium dynam- 
ics during the condensation of a spin-1 Bose gas with 
SOC. In particular, we focus on the formation of topo- 
logical defects in the limit of rapid temperature quench. 
According to Kibble-Zurek mechanism [13 [18], topolog- 
ical defects can be created via phase transitions at finite 
temperatures, which are caused by spontaneous sym- 
metry breaking and thermal fluctuations near the crit- 
ical point. By solving the stochastic projected Gross- 
Pitaevskii equation (SPGPE), we show that, in a spin-1 
BEC, the combination of SOC, spin-exchange interaction 
and thermal fluctuations can generate meron-like excita- 



tions. A meron is a peculiar topological defect that was 
originally hypothesized as a half-instanton in the parti- 
cle physics [19] and later a half-skyrmion in the quan- 
tum Hall systems [20j [21] . Since merons carry half unit 
of topological charge, it is believed that isolated merons 
can only be observed when particular boundary condi- 
tions are imposed. By far, merons have been created 
in the superfluid 3 He- A [22j [23] in a rotating cylinder, 
and the spinor BEC with a special constraining magnetic 
field [24 . In what follows, we show that, stable collec- 
tive excitations such as the crystalline orders of merons 
and other isolated variants can be created in a rapidly 
quenched spinor BEC with SOC. 

The order parameter of a spin-1 BEC is given by 
^ = (^i, ^o, ^-i) T , where T stands for the trans- 
pose and ^ mF (rriF = ±1,0) denotes the macroscopic 
wave function of the atoms condensed in the spin state 
llj/np). The total particle number, A 7 ", and total mag- 
netization, M, are normalized by J \^\ 2 d 3 r = A/", and 
/ (j^i| 2 - |^-i | 2 ) d 3 r = M. In the following, we con- 
sider SOC of the form H so — ^ a V a p a F a , where V aj p a 
and F a are respectively the coupling strength, the linear 
momentum and the 3x3 matrix of the spin-1 angular 
momentum in the a (= x, z) direction. In the absence 
of magnetic field, the dynamics of \l/ is described by the 
following coupled nonlinear Schrodinger equations, 
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Here H = -Ti 2 V 2 /2ra + U(r)+ g n |\l>| 2 denotes the spin- 
independent part of the Hamiltonian, with U(r) being 
the trapping potential. The coupling constants g n and g s 
characterizing the density-density and spin-exchange in- 
teractions, respectively, are related to the s-wave scatter- 
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ing lengths ao and CL2 in the total spin channels F tota i = 0, 
2 by g n = Anh 2 (ao + 2*22) /3m, g s = 47ra 2 (a2 — ao) /3m 
[25] [26]. Note that g n > 0, whereas g s can be either pos- 
itive or negative. A spin-1 BEC is said to be ferromag- 
netic (FM) when g s < 0, and antiferromagnetic (AFM) 
when g s > 0. As we shall focus on the dynamics of 
spin texture, we introduce the basis set ty a (a = x, y, z), 
such that \P±i = (±V X +i i &y) /V% and ^ = ^z- As 
a result, F a \a) = 0, and the spin texture, which is 
parallel to the local magnetic moment, is defined by 
S (r) = i& x */ |^| 2 where ^ = (V x , * y , * z f [27 . For 
later use, we define the unit vector s (r) = S (r) / |S (r)|. 
In the mean-field theory approximation, the dynamics 
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FIG. 1: (color online) Snapshots of (a) |^-i| , (b) |^ |, 
(c) |^i | 2 , and (d) |*| 2 of an 87 Rb (FM) spinor BEC during 
quench with V x = 1, V y = 2. The rightmost column shows 
the phase profile of equilibrium state. The particle numbers 
in the equilibrium state are N±i « 3.65 x 10 3 , iVo = 7.32 x 10 3 . 

of a BEC at nonzero temperatures can be described by 
the SPGPE [28-30 , based on the assumption that the 
system can be treated as a condensate band in contact 
with a thermal reservoir comprising of all non-condensed 
particles. The condensate band is described by the trun- 
cated Wigner method [31] including the projected c-field 
method, while the non-condensate band is by the quan- 
tum kinetic theory [32] [33]. A direct generalization to 
the spinor BECs with SOC leads to the following set of 
coupled SPGPEs [34] 
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FIG. 2: Snapshots of the vector field s (r). Arrays of merons 
of the same handedness can be seen in (b). All merons drift 
away towards the periphery of the condensate as shown in (c) 
and (d). 



where T and /i denote the final temperature and chem- 
ical potential, jj the growth rate for the j-th compo- 
nent, and dWj/dt is the complex- valued white noise 
associated with the condensate growth. The projec- 
tion operator V restricts the dynamics of the spinor 
BEC in the lower energy region below the cutoff energy 
En. For simplicity, we shall consider a 2-dimensional 
isotropic trap, U(r) = muo 2 (x 2 + y 2 ) /2. The numeri- 
cal procedures for solving SPGPEs, are described as fol- 
lows. First, the initial state of each Sfrj is sampled by 
using the grand-canonical ensemble for free ideal Bose 
gas at a temperature To below the critical tempera- 
ture and of chemical potentials /i- j0 . The spatial de- 
pendence of the initial state can be specified as a lin- 
ear combination of plane waves with discretized mo- 
mentum k = 27r (n x ,n y ) /L (n x , n y G Z and L is the 
size of the system), i.e., ^ (t = 0) = Xlk^ ^k^k ( r )> 



where ^ k (r) = e* k r . The condensate band lies below 
the energy cutoff Er > E\^ - 
but ion is sampled by a^k = (^Y?,k 
N jM = [exp((£ j?k - fjL jm0 ) /k B T ) 



h 2 |k| 2 /2ra. The distri- 
1/2) 77 - jk where 

l] 1 and 77 jk are 
the complex Gaussian random variables with moments 

(Vj,\cVj,v) = {nl\tVj,w) = and (vj^Vlv) = <W- 
Second, to simulate the thermal quench, the tempera- 
ture and chemical potential of the non-condensate band 
are altered to the new values T < Tq and \i > fij . For 
convenience, we adopt the oscillator units in the numer- 
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ical computations, and the lengt h, time and energy are 
respectively scaled in units of ^Ti/muo, cj _1 and Tluj. 
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FIG. 4: (Color online) (a) The 3-dimensional orientation of 
the spin textures in Fig. 1. The color indicate the magnitude 
of s z . (b) The topological charge density a of the equilibrium 
spin texture, (c) The corresponding spatial distribution of 
|S(r)|. 



FIG. 3: (color online) Snapshots of (a) |^-i| 2 , (b) |^o| 2 , 
(c) |^i| 2 , and (d) of the 87 R b (™) s P in or BEC at zero 
magnetization with V — 1.8. The rightmost column shows the 
phase profile of equilibrium state. The corresponding particle 
numbers are N ±1 « 3.56 x 10 3 , iV = 7.08 x 10 3 . 

We first study the condensation of 87 Rb (g s < 0) with 
SOC. The initial states are sampled at /c^Tq = 500, 
fi j = 2. Here we set E R = 50.5, k B T = 0.067, /i 25, 
and h^j/kBT = 0.05. In this paper, we consider the in- 
plane coupling, i.e., V x ,V y ^ and V z = 0. We begin 
with the cases of \V X \ ^ \V y \. In Fig. 1, stripe structure 
develops in the density profile of each component dur- 
ing condensation, but fade away gradually. On the other 
hand, the equilibrium phase profiles become periodically 
striped along one direction, which is is exactly the plane 
wave (PW) state found in Ref. [11 . The evolution of the 
spin texture is shown in Fig. 2. Clearly, spin domains 
form shortly after the quench starts. The domain walls 
comprise of sheets of spin vortex which drift outward 
and disappear while the system approaches equilibrium 
marked by a uniform spin alignment on the xy plane. 
The PW state is bound to occur when \V X \ ^ |V^|, or 
\V X \ = \V y \ = V < 0.8. When \V X \ = \V y \ = V > 0.8, 
periodic structures, which are caused by the formation of 
grids of dark soliton in \£±i and vortex lattice in \I>o, may 
appear in the equilibrium density profiles of all compo- 
nents in Fig. 3. In Fig. 4(a)-4(b), the spin texture con- 
sists of two interlacing square lattices of spin vortex with 
opposite vorticity, and we denote it as the spin vortex 
lattice (SVL) state. When V > 0.8, PW and SVL state 
appear alternatively with definite probabilities depend- 



ing on the magnitude of V. Using the imaginary time 
propagation method, the SVL state is shown to have a 
higher energy than the PW state. In the SVL state, all 
spin vortices have similar structure — the central spin 
always orients to the z axis, while the others increasingly 
twist and finally lie on the xy plane, forming a circulation 
pattern away from the center. The topological charge 
density, a = s • (d x s x d y s) /4-7T, is plotted in Fig. 4(b) 
and the vortex with left /right-handed circulation has a 
positive/negative a. This is just the Mermin-Ho vortex, 
or meron [35 j. Integrating Q = J crdPr over a primitive 
unit cell, we identify Q = ±1/2, corresponding to merons 
and antimerons, respectively. Due to the FM nature of 
the condensate, a meron and an antimeron will pair up to 
form a vortex dipole, which was predicted in the bilayer 
quantum Hall system as the lowest energy excitation [2T] 
and in a fast-rotating highly spin-polarized spinor BEC 
[27] and now acts as the building block of the SVL. For 
conciseness, in the following, we shall not distinguish the 
merons from antimerons unless specially noted. 

To gain more insight into the SVL state, we plot the 
equilibrium momentum wave functions, (p), and find 
that they are all sharply peaked at p ~ ±2.8e x , ±2.8^. 
This suggests that the equilibrium state might be con- 
sisted of 4 plane waves with p = ±q, ±q ; , where q q' = 
and |q| = |q'|. The counter-propagating modes with 
p = ±q, ±q' form two orthogonal standing waves that 
superimpose to generate the periodic structures in all 
density profiles in Fig. 3. In the absence of trapping 
potential, it was shown that, the one-particle Hamilto- 
nian, p 2 /2m + Vp • F, is minimized by |p| =V [HI [12] . 
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and hence 
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where = tan 1 (q y /q x ). Consequently, the 
are given by, 



spin textures 



S (r) = (cos cos u sin v + sin sin u cos 

cos sin it cos v — sin # cos u sin cos u cos v) (5) 

where u = V(xcosO + ysinO) and v = V(—xsm9 + 
2/ cos 0) . The above analytical results fairly reproduce the 
spin texture plotted in Fig. 4(a) with an effective V 
in the presence of trapping potential. From Fig. 4(c), 
we see that the spin density vanishes at the center of a 
vortex quadruple (4 mutually adjoining merons), which 
happens to be the density troughs of that is filled 
by the particles of \I/o- The meron and polar cores are 
centered at r meron / polar =# ((9) d meron/polar where R (0) is 
the rotation matrix on xy plane, d meron = (n,/)7r/V, 
dpoiar = (n + 1/2, 1 + 1/2) tt/V, and n, I e Z. It should 
be noted that the cores of the meron are located at the 
vortices of \P . With Eqs. ([U-Q, the energies of PW and 
SVL states for a homogeneous spinor BEC can be calcu- 
lated analytically. As PW and SVL states both yield the 
same minimized energy for the single particle Hamilto- 
nian with a fixed TV, we only need to consider the inter- 
acting energies. It is straightforward to show that both 
PW and SVL states lead to the identical density-density 
interaction energies and thus the only difference is the 
spin-exchange term. As a result, the spin-exchange in- 
teraction energies per unit area, e sp i n = g s N j S 2 (r) <i 2 r, 
are g s N and 3g s N/4 for PW and SVL states, respectively. 
In the current case, g s < 0, and thus the SVL state has 
a higher energy. 

The SVL state is sustained by the vorticity originated 
mostly from SOC. As \i rather than N is fixed in our 
simulations, increasing V will increase N and the mass 
vortices in SPq. Since PW and SVL states are gapped 
by an amount of energy proportion to A 7 ', it is expected 
that there exists a threshold V c , beyond which the SVL 
will barely appear. This implies that the PW state to- 
tally prevails in the large V limit. The value of V c can 
be estimated by considering a homogeneous spinor BEC, 
in which the core size of the mass vortex nucleated in 
the vfip = component is £ ~ fi>/ |^o| V2mg n . The in- 
stability of SVL state occurs when the lattice constant 
is comparable to £ , such that a polar core will par- 
tially overlap with the cores of contiguous merons, giv- 
ing V c ~ 7T 

|^o | y/2mg n /h. By the same token, if the size 



of the condensate is smaller than the lattice constant, 
the condensate does not accommodate SVL, which oc- 
curs when V < 0.8 in our previous result in a trap. In 
the limit of vanishing SOC, the spatial periodicity in the 
phase profiles of PW and SVL states is infinitely pro- 
longed, leading to uniform distribution in all phase pro- 
files — the manifestation of a normal spinor BEC. We 
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FIG. 5: (color online) Spin textures and the associated distri- 
bution of \S (r)| in a 23 Na (AFM) spinor BEC for the inverted 
merons with (a), (b) M/N = 0.88; (c), (d) M/N = 0.78. The 
particle numbers in the equilibrium state are Ni ~ 4.66 x 10 4 , 
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5.14 x 10 3 , N-! & 3.74 x 10 2 in (a) and Ni n 4.18 x 10 4 
9.51 x 10 3 , N-i « 9.03 x 10 2 . 
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FIG. 6: (Color online) (a) The stereographic projection of the 
spin textures onto S 2 for (a) a meron; (b) an inverted meron. 

now consider the spinor BEC of 23 Na (g s > 0) with SOC. 
Using the same parameter setting, we conclude that if the 
spin-1 gas is unpolarized initially, no spin texture will be 
produced in the condensate, a conclusion already given 
in Ref. [11 . The equilibrium spin texture for a spin- 
polarized state with M/N & 0.88 is shown in Fig. 5(a), 
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where all spins in the outer region point to the +z di- 
rection, tilting gradually, and finally lie on the xy plane 
while approaching the center region. Consequently, a 
polar core is formed at the center, with \I/o filling in the 
vortices of both \£±i which are separated but locate prox- 
imately to the center. Such a configuration corresponds 
to a charge of Q = — 1/2. In Fig. 6 the stereographic pro- 
jections for a meron and the texture of Fig. 5(a) are plot- 
ted. Clearly, Fig. 6(b) is exactly the inverted figure of 
Fig. 6(a), and thus we term the structure of Fig. 5(a) the 
inverted meron. In Fig. 5(c)-5(d), three inverted merons 
with distorted cores are formed with M/N ^0.78. It is 
easy to verify that the single isolated inverted meron has 
a lower energy. 

In summary, we have investigated the non-equilibrium 
dynamics of spin-1 BECs with SOC in the limit of rapid 
quench. Crystallization of merons and polar core vortices 
are predicted to arise in the FM spinor BEC. Likewise, 
isolated inverted merons can be created in the highly po- 



larized AFM spinor BEC. Our studies provide a method 
to create nontrivial structure of merons and thus an op- 
portunity to probe into the fundamental properties of 
meron-like matter. Following the experimental meth- 
ods in Ref.[2], our predictions can be realized in prin- 
ciple, except in lifting the degeneracy of the hyperfine 
spin states of F = 1, a weaker magnetic field is needed 
to avoid the decoupling of spin states due to quadratic 
Zeeman shift. The experiment may start by trapping 
a thermal spin-1 Bose gas and then rapidly quench it 
to the quantum degenerate regime. Finally, the equilib- 
rium spin texture can be resolved in situ by using the 
polarization-dependent phase-contrast technique [36] . 
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